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The discovery of discrete sources of x-rays'^ in the 

sky has led to much speculation as to the mechanisms 

2 — 6 

responsible for the x-ray emission . Among these 

suggestions has been the thermal emission of x-rays from 

the hot surfaces of neutron star remnants of supernova 
2 3 

explosions ' . Calculations of the cooling of such a 

neutron star following its formation have seemed to 

indicate that the temperature would still be several 

million degrees at an age of several thousand years. 

Such cooling calculations have been based upon estimates 

of photon emission from the surface and of neutrino- 

antineutrino pair emission by the plasma process from 

the interior. There have recently been new suggestions 

regarding other possible mechanisms of neutrino emission 

7 

from the interior , and this has led to the impression 

that neutron stars might cool too rapidly to be observa- 

5 8 

ble as celestial x-ray objects ' . We believe this view 
to be overly pessimistic, and we give here a brief summary 
of some of our cooling calculations for neutron stars. 

The neutron star models we have used in our calculations 
have been based upon composite equations of state with 
nuclear forces included. The details of the construction 
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of the equations of state have been described in a thesis , 
and will be published in due course. At the lowest 
densities, but still at substantially high temperatures, 
the composition was assumed to be iron. As the density 
increases, the main contribution to the pressure comes 
from the degenerate electrons, and the mean molecular 
weight per electron gradually increases as the high Fermi 
level of the electrons forces the nuclear composition of 
the material to change toward higher mass number. At 
densities between 10^ and 10^ gm/cm^ the heavy nuclei 
gradually dissolve into neutrons. The system then 

consists of degenerate neutrons, protons and electrons. 

15 3 

Near and above 10 gm/cm , many other particles appear 

in the mixture, commencing with the p meson 

and the £ hyperon. In constructing our equations of 

state, the degeneracy pressures of all the individual 

fermions were added, and nuclear force terms were also 

included. 

Nuclear forces are generally attractive at large 
internucleon distances and repulsive at small internucleon 
distances. In this work we chose forms of the potential 
interactions between neutrons suggested by Levinger and 
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Simmons . Their potentials V and V were utilized. 

B Y 

Both of these potentials are attractive at low densities, 

although V is somewhat more attractive V . At hiqh 
Y ' B 

densities rapidly turns repulsive, while only 
slowly turns repulsive. For our neutron star equations 
of state, we have applied these potentials between the 
baryons without distinction as to the type of baryon. 

At densities less than or equal to nuclear densities, 
the character of nuclear forces is reasonably well known 
and is given to a rough approximation by either of these 
potentials, and the composition of the matter is mostly 
neutrons, for which the potentials were originally 
constructed. At much greater than ordinary nuclear 
densities many different types of baryon are present, 
and the rapidity with which nuclear forces turn repulsive 
is very speculative. Hence the two potentials and 
tend to span a range of possible behavior of the 
nuclear forces at high densities, and the differences 
in the neutron star models which result from the adoption 
of one or the other potential will give an indication of 
the uncertainty due to lac|£ of knowledge in this area of 


physics . 
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In the case of a perfect fluid there is a general 
relativistic limitation on the pressure such that it 
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cannot exceed one-third of the proper energy density 
In a fluid with suitable anisotropic properties, this 
condition may be violated, but the ultimate relativistic 
condition still remains that the pressure cannot exceed 
the proper energy density'*' 1 . If this were to be violated 
the speed of sound would exceed the speed of light in the 
medium. Accordingly, the composite equation of state 
constructed as described above was cut off with one of 
these two pressure saturation conditions. 

The general relativistic equations of hydrostatic 

12 

equilibrium were derived by Oppenheimer and Volkoff 
For a given assumption about central density, the numerical 
integration of the differential equations of hydrostatic 
equilibrium gives models with uniquely determined masses. 
The gravitational and proper masses determined in this 
way as a function of central density for the two composite 
equations of state are shown in Figure 1. Several 
interesting comments follow from this figure. It is 
evident that the details of the structure of the neutron 
star models are very sensitive to the uncertainties in 
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the rate at which nuclear repulsive forces enter in the 
baryon mixture at high densities. For each equation of 
state the mass rises with increasing central density 
toward a principal peak, beyond which it falls and then 
oscillates. It is interesting to note that the gravita- 
tional mass is less than the proper mass over the entire 
range of neutron star models, thus showing that such 
models are gravitationally bound. This is contrary to 
the behavior of models constructed with non-interacting 
neutron gases. An expanded discussion of the stability 

of the models beyond the main peak will be given 
13 14 

separately ' . It may also be noted that pressure 

saturation does not set in until the vicinity of the 
principal peak has been reached. Hence uncertainties 
in the pressure saturation effects play no significant 
role in the discussion which follows. 

In order to calculate representative cooling curves 
for neutron star models, three models were chosen so 
that one was of low mass at the low density base of the 
principal peak, one was halfway up the peak, and the 
other was near the top of the peak. To these models 
were fitted hot atmospheres corresponding to a series 
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of surface temperatures . These atmospheres were composed 
of both iron and magnesium; the results were nearly the 
same and only the neutron stars with iron envelopes are 
discussed here. It was typically found that the temperature 
rose by about a factor 50 between the photosphere of the 
neutron star and the interior where the high thermal 
conductivity of the electrons assured a flat temperature 
distribution. The atmospheres were constructed with the 
help of opacities calculated from the Los Alamos opacity 
code of A.N. Cox and his collaborators. The atmospheric 
structure was determined by requiring that the luminosity 
should not change from one layer to the next. 

The heat capacity of the neutron star models is a 

15 

function of their temperature . The presence of 
nuclear forces in the equation of state will modify the 
heat capacity by an amount which typically can be of the 
order of a factor 2, according to rough estimates which 
we have made. We did not take such modifications into 
account in making the actual cooling calculations. 

The cooling of the neutron stars is due to the 
combination of neutrino emission from the interior 
and photon emission from the surface. The construction 



of the envelope automatically provided us with the 
photon cooling rates. Three neutrino cooling rates 
were taken into account, as follows: 

(1) Neutrino pair emission from the plasma process. 

These neutrinos arise from the decay of plasmons in the 

degenerate electron gas in the interior of the neutron 

star. The rates have been given by Adams, Ruderman, 

16 

and Woo and by Inman and Ruderman 

(2) The URCA process. Bahcall and Wolf have 

recently given an estimate of neutrino emission from 
0 

this process , in which neutrons decay into protons 
and protons capture electrons. The rate is somewhat 
greater than that of the plasma process. 

(3) The neutrino bremsstrahlung process. The 
neutrino pairs are emitted when electrons scatter from 
positive or negative baryons in the interior of the 
neutron star. Ruderman and Festa have kindly provided 
us with the following approximate preliminary expression 


for this process : 
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where Z is the effective charge of the electron 
scattering centers, n is the number density of such 
centers, and n we take here to be the bar yon number 
density. Ruderman and Festa have suggested that there 
may be proton clustering in neutron star interiors in 
the presence of very large numbers of neutrons, so that 
the effective charge of scattering center might be 2. 
However, it may well be that under conditions of interest 
in the interiors of neutron stars there will also be 
large number of E hyperons, which may hinder the 
clustering process. Consequently, we have chosen to 
take the effective charge of a scattering center equal 
to unity and to count as the scattering centers both 
the protons and the E hyperons. This process is less 
important than the URCA process at high temperatures, 
but it is more important at low temperatures . 

The cooling curves for the six chosen neutron star 
models as a function of age are shown in Figure 2. It 
may be seen that the rate of cooling has a significant 
dependence on the mass of the star. The low mass stars 
cool quite rapidly, but the medium and heavy mass stars 
still have temperatures exceeding 2xl0 6 °K for times of 
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the order of 10 or 10 years . Hence it is evident 
that thermal emission of x-rays from neutron star 
surfaces should continue to be regarded as candidates 
for identification with some of the x-ray sources which 
are being found in the sky. 

Bahcall and Wolf have raised the question of neutrino 
cooling from pion decays in neutron star interiors. Such 
pion decays can occur only if pions should have a small 
effective mass in the presence of a largely neutron gas. 
Both Bahcall and Ruderman have indicated to us (private 
communications) their expectations that, under the 
conditions in which pions may be present in a neutron 
star, there will be a predominantly repulsive interaction 
between the pions and the neutrons . This would raise 
rather than lower the effective mass of the pions, and 
make it very unlikely that pions will be present in the 
interiors of neutron stars on the low density side of 
the principal peak. 

The remarkable NRL rocket experiment carried out 
during a lunar occultation of the Crab Nebula has shown 
that the x-ray source associated with the Crab Nebula 
has dimensions very much larger than those of a neutron 
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star. One of us has suggested that these x-rays may be 
due to the synchrotron process from high energy electrons 
accelerated in the magnetosphere of a vibrating neutron 

g 

star . Similar nonthermal processes may well be 
associated with other x-ray sources, if they are neutron 
stars, and hence nonthermal components to x-ray spectra 
may be common. Hence we believe that many more highly 
refined experiments will be necessary before the true 
nature of the celestial x-ray sources will be determined. 

We are grateful to Dr. Ruderman and Mr. Festa for 
communicating to us in advance of publication their 
preliminary results on neutrino emission by the 
bremsstrahlung process. 
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FIGURE CAPTIONS 


Figure Is Gravitational and proper masses of the neutron 
star models constructed with both composite equations 
of state and plotted as a function of central density. 
The different models resulting from the two pressure 
saturation conditions at the highest densities are 
separately indicated. 

Figure 2: The cooling curves for selected neutron star 

models fitted with atmospheres of iron, with three 
models for each composite equation of state. In the 
portions of the curves to the left of the crosses 
cooling occurs predominantly by neutrino emission 
from the interior? to the right of the crosses photon 
cooling from the surface predominates. 
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